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Using the Landshoff-Nachtmann two-gluon exchange model of the pomeron, 
the double pomeron exchange contribution to production of gluon pairs in 
the central region of rapidity is calculated. The results are compared with 
those for production of quark- ant iquark pairs. 

1 Introduction 

Production of the Higgs boson by double pomeron exchange was recently studied by 
several authors [1]. Unfortunately the published theoretical estimates of the cross- 
sections are widely different and thus do not permit to obtain reliable predictions 
needed for future experiments. This reflects our present limited understanding of the 
nature of the diffractive (pomeron) reactions. 

One way to reduce this ambiguity is to calculate the cross-section for other double 
pomeron exchange processes and compare them with existing data. In the present 
paper we study this problem using the Landshoff-Nachtmann model of the pomeron [2] 
to calculate the production of two large transverse momentum gluon jets shown in Fig. 
1. In this way we hope to obtain a better estimate of the Higgs production studied in 
this model some time ago [3]. 

Our calculations follow closely the method used in [4] [5] where the cross-section 
for production of heavy quark-antiquark pairs was calculated. The results for two 
gluon jets, together with those for quark jets, allow to obtain the full cross-section for 
double-diffractive jet production in the Landshoff-Nachtmann model. 



1 e-mail: bzdak@th.if.uj.edu.pl 

1 




Figure 1: Production of two gluons by double pomeron exchange. 

2 The relevant diagrams 

In the Landshoff-Nachtmann model [2] the pomeron is described as an exchange of 
two " nonperturbative" gluons, which takes place between a pair of quarks of colliding 
hadrons. Thus to calculate the cross-section for reaction pp — ► pp + gg by double 
pomeron exchange we have to take into account twelve diagrams presented in Fig. 2. 




Figure 2: The diagrams which have to be taken into account in the process of the 
gluon pair production by double pomeron exchange. The dashed lines represent the 
exchange of the nonperturbative gluons. 

Other diagrams, like these presented in Fig. 3, either do not contribute significantly 
because of large transverse momentum on one of the quark internal lines or simply 
vanish because pomeron is the colour singlet. 
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Figure 3: The examples of the diagrams which do not contribute. 

Following [3] we assume that up to highest powers of s/s\ and s/s2, the total 
amplitude given by diagrams in Fig. 2 is of the form: 

Here a± = a (ti) , a 2 = a (t 2 ) and a (t) = 1 + e + a't is the pomeron Regge trajectory 
with e w 0.08, a' = 0.25 GeV~ 2 ( si, s 2 , ti, t 2 are marked in Fig. 1). So, in order to 
calculate the amplitude for production of two gluons by double pomeron exchange, it 
is enough to take diagrams from Fig. 2 at ti — t 2 — 0. 

Even with this simplification, however, the direct evaluation of all diagrams of Fig. 
2 is still tedious. Fortunately, as shown in [6] and exploited in [3], it turns out that 
up to terms of the order e, the matrix element Mfi\ tl=t2=0 is equal to the sum of the 
following three diagrams: 



p, q, k 1 p, q 1 k. Pi q, k 1 




Figure 4: Putting the quarks lines q\ and q 2 on shell is equivalent to calculating the 
amplitude for production of gluon pair by double pomeron exchange at t\ = t 2 = 0. 

where the lines q\ and q 2 are put on shell and t 1 — t 2 — 0. 

The calculation is performed in two steps. First, we calculate Mfi\ tl= t 2 =o for col- 
liding quarks, given by the sum of the diagrams presented in Fig. 4 and use (1) to 
introduce the Regge behavior. To obtain the cross-section the presence of three quarks 
in the proton must be taken into account. 

3 Matrix element 

We take the masses of the light quarks to be approximately zero and, following [3], use 
the Sudakov parametrization for momenta shown in Fig. 4: 
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x y 

Q = ~Pl + ~P2 + V, 

s s 

h = Xipi H P2 + U1, 

s 

&2 = — Pi + 2/2P2 + V 2 , 
S 

r 2 = x g px + y g p 2 + v g , (2) 

where v, vi, v 2 , v g are two-dimensional four-vectors. 
Let us express Mfi\ tl=t2=0 as the following sum: 

M fi \ tl=t2=0 = Mjf + Mjt\ (3) 
where is the sum of the first two diagrams presented in Fig. 4, and M^f being last 



diagram in Fig. 4. Following closely the method used in [5] we obtain (see appendix): 

Mf = -A^/j^nrse^, (4) 

M$ = -A^5%5 UW J d 2 vf (v) exp (=p) , (5) 



where: 



f ^ = ~2 ( s2e i £ 2 (S2 - y g ) %g ~ 2se 1 ve 2 v + 2se 1 ve 2 p 2 (5 2 - y g ) - 2.SE 2 VEip 2 (5 2 - y g )) 
I 

+ -z^ (s 2 £iE 2 (5 1 - x g ) y g - 2se 1 ve 2 v + 2se 2 ve 1 p 2 y g - 2se 1 ve 2 p 2 y g ) . (6) 

Here A = i(D G 2 ) 3 g 2 /(27i) 2 G 2 , G and g are the non-perturbative and perturbative 
quark-gluon couplings, r = fi 2 / (1 + X\ + y 2 ). D (p 2 ) = D Q exp (p 2 / fi 2 ) with /j, m 1 
GeV and G 2 Dq m 30 GeV _1 /x _1 , is the non-perturbative gluon propagator. By 5\ and 
5 2 we denote 1 — x\ and 1 — y 2 respectively. The indices u and w denote the colour of 
the produced gluons. S\ and e 2 are the polarization vectors of gluons chosen to satisfy 
the relations: 

^er(A)<(A) = -^ + Ar + Krf i = l,2. (7) 

A=l,2 PlTi 

So that, in addition to the constraints E\T\ = e 2 r 2 = 0, we have further constraints: 
£iPi = e 2 px = 0. 

Performing the integration in (5) is a rather difficult task. Following [5] we expand 
the integral over d 2 v in power series: 

//— V 2 \ ITT 2 
d 2 vf (v) exp ( — — J =anr + Y^ ~ a n + •••> ( 8 ) 
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where a and Qjjj £11*6 coefficients of expansion of / (v): 

ld 2 f (v) 

a = f(v = 0), %, = 2^W k= °- (9) 

One can check that in the first term of (8) and M^f exactly cancel each other. 

So we are forced to calculate the second term. After rather lengthy calculations we 
obtain: 

4 f 2 f 1 2 5 2 I 

M fi tl=t2= o = -A— ^y uw 7rr <^ eie 2 H -^iPa^ f • (10) 

12 t x aVg sx g y g 6 l J 

Taking the square of the amplitude (10), averaging and summing over spins, colours 
and polarizations one arrives at the simple formula: 



\M filtl=t2=0 \ - — — 



2 = # 



^2\3 ,,4\ 2 / „2 \ 2 



where: 



(12) 2 ^ 9(2tt) 2 y VC 2 
Factors 8 and (y^) 2 are related to the summing over colours of produced gluons and 

(3) 

to the square of colour factor of diagrams represented by , respectively. 

Taking into account (1) where (s/si) = (l/5 2 ), (s/«2) = (l/^i), £i = — # 2 and 
t 2 = — v\ we obtain: 



= 7 , 2 , AA vfe^c**, ex P W (*i +*»))■ (12) 

(<W 5 i 5 2 

The factor exp (2/3 (ti + t 2 )) takes into account the effect of the momentum transfer 
dependence of the non-perturbative gluon propagator with (5 — 1 GeV -2 [7] [8]. 



4 Cross-section 



2 

Having calculated |M/j| we can write the formula for the cross-section: 



81 



2s (2tt) 8 2! 



y (M^fF^F^)] 2 ^//, (13) 



where the factor 81 takes into account the presence of three quarks in each proton, 2! is 
an identical final state particle phase space factor 2 and F (t) is the nucleon formfactor 
approximated by: 

F (t) = exp (At) (14) 



2 We would like to thank J. R. Cudell for a correspondence about this point. 



5 



with A = 2 GeV 2 . Differential phase-space factor dPH has the form: 



dPH = d 4 h5 (k\) d 4 k 2 S (k 2 2 ) d A ri 5 (rf) d 4 r 2 5 (r 2 2 ) x 

(fc°) (kg) (r°) (r 2 °) 5^ ( Pl +p2 -k 1 -k 2 -r 1 - r 2 ) . (15) 

It turns out that a substantial part of the integrations can be performed analytically 
[5] . Denoting by -E^in the minimum value of transverse energy of the produced gluon 
(we integrate over E J > E^ in ^ 0) and by A the maximum value of the energy loss of 
the initial hadrons (Sip < A ~ 0.1), we obtain: 

MA) 

2x 

r2 



a (ETJ =cjdxx(l- x 2 f (in i±£ + ^ 
o 

In »' ? , (16) 



where /i (A) = ^1 - 5 2 /A 2 , 5 2 = 4(ElJ 2 /s and: 



C = 9tt 



7T 2 (DoG 2 ) V 4 \ f<? 



ll8(2vr) 6 ^ in a'/ V^ 2 



.2 \ 2 



The major uncertainty in above result is the value of the non-perturbative coupling 
constant G. 

Finally, let us remind the result for quark- ant iquark production by double pomeron 
exchange [4]: 



h(A) 

a UkJ niQ ) 2 ) = C [ dxx ( 1 - m 2 }~* - ) (1 - x 2 Y +2e ( ln — + -^-] 
Ky mmJ ' J V rn 2 + (Knn) 2 J ' \ l-x l-x 2 ) 



1 ^ + ln^!l 



^ - vzg , ; 5 , (17) 



+ I In i^i h±l + i n J, ' 



where ft, (A) = ^/ 1 — 5 2 /A 2 , 5 2 = 4 (m 2 + (fc^J 2 ) /s, m is the mass of the produced 
quark, and: 



~ 1 / s \ 2t n 2 (D G 2 y^m 

C = — 7T 



2 ,2X2 



3" V4(m2 + (A£JV V 18 ( 27r ) 6 «'( m2 + ( A; m i „) 2 )y VG 2 , 
By /c^ in we denote the minimum value of transverse momentum of the produced quark. 

5 Numerical results and discussion 

Let us remind the following numbers that appear in (16) and (17): e = 0.08, a' = 0.25 
GeV" 2 , fi = 1 GeV, G 2 D = 30 GeV-V 1 , A = 2 GeV" 2 , (3 = 1 GeV~ 2 , A = 0.1. 
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The minimum value of transverse energy of the produced gluon E^ in is chosen at 5 , 7 
and 10 GeV . At Tevatron energy, y/s = 1.8 TeV, the results for the cross-section a 
(multiplied by (G 2 /4tt) ) are as follows: 



[TeV] 


= 5 GeV 


£T m = 7 GeV 


Ej niQ = 10 GeV 


1.8 


(G 2 /47r)" a = 3 fib 


(G 2 /4tt)V = 0.9 fib 


(G 2 /47r)" a = 0.3 ^b 



Table 1: The numerical results for the cross-section for production of gluon pairs by 
double pomeron exchange. The calculation is performed at Tevatron energy y/s = 1.8 
TeV. 

The running coupling constant g 2 /47r was evaluated at 2E^ in , i.e. 0.17, 0.15 and 
0.14 for El in = 5, 7, 10 GeV respectively. 

The dependence of the cross-section versus the center mass energy Eqm — \/s for 
the range of y/s taken from 0.5 TeV to 40 TeV is shown in Fig. 5. 




Figure 5: Double pomeron exchange contribution to the cross-section for the gluon 
pairs production process versus the center mass energy. 

In order to obtain the total cross-section for production of jets by double pomeron 
exchange, we have to consider the quark- ant iquark production process. It is obvious 
that for the case in which the minimum value of transverse momentum of the produced 
quark k^ oin = 0, the relations hold: a uU ~ <r d j > a sS » a C c > &bi> > a tt- The values 
for o"cc, cr b i and a t i are determined by (17) (this expression is reliable for heavy quarks 
production process). For more details see [4] [5]. 

In the experiment we observe jets with some minimal E~^ ni ~ GeV, however, and 
these relations are not valid. At Tevatron energy y/s = 1.8 TeV and E^ = 7 GeV we 
have: 
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[TeV] 


(G 2 /4tt)V cc - [nb] 


(G 2 /47r)V 5 [nb] 


(G 2 /4n) z a tI [nb] 


1.8 


0.3 


3.1 






Table 2: The numerical results for the cross-section for production of heavy quark 
pairs by double pomeron exchange. The calculation is performed at Tevatron energy 
y/s = 1.8 TeV and = 7 GeV. 

As the masses of the heavy quarks we take m c = 1.2 GeV, m b = 4.2 GeV, m t = 175 
GeV (there is not enough energy to produce a tt pair). The running coupling constant 
g 2 /4:7r was evaluated at 2E£ lin , i.e. 0.15, 0.15 and 0.1 for cc, bb and tt respectively. 

The dependence of the cross-section on the center of mass energy Eqm — \fs for 
the range of y/s taken from 0.5 TeV to 40 TeV is shown in Fig. 6. 




Figure 6: Double pomeron exchange contribution to the cross-section for production 
of heavy quark pairs process versus the center mass energy. 

As we see a b i > a c5 . This is a consequence of the fact that we observe the jets 
with minimal E^ = 7 GeV . Since there are very few light quarks with E^ in ~ GeV, 
the cross-section for production uu dd ss pairs with such minimal energy is negligible. 
To summarize, the cross-section for quark-antiquark production with minimal Ej nin ~ 
GeV is fully dominated by the cc and bb pairs production. 

One sees that qq contribution to production of jets is about three orders of magni- 
tude smaller than that of gg contribution. Similar observation was first made in Ref. 
[9]. 

Our cross sections are calculated for the case where both initial protons are scattered 
quasi-elastically. If one allows for diffractive dissociation of the incident hadrons i. e. if 
one puts A = in (14) the cross sections increase about half order of magnitude. 



8 



6 Conclusions 



We have calculated the cross-section for gluon pair production by double pomeron 
exchange in the Landshoff-Nachtmann model. The obtained results were compared 
with those for production of quark- ant iquark pairs calculated in the same model. It 
turns out that cross-section for jet production is strongly dominated by gg production. 

I would like to thank Professor Andrzej Bialas for suggesting this investigation and 
helpful discussions. Discussions with Prof. Maciej A. Nowak, Prof. Jacek Wosiek, 
Prof. Michal Praszalowicz and dr Leszek Motyka are also highly appreciated. 

7 Appendix 

In the present appendix some main steps leading to the expressions (4) and (5) are 
more explicitly shown. 

Let us denote the colour structure of the diagrams represented by ilfif as follows: 



d e f d e f 




9 h k 9 h 



Figure 7: The capital letters denote colour indicies of the non-perturbative gluons and 
small letters colour indicies of quarks and perturbative gluons. 

Performing colour calculations for the first diagram in Fig. 7 (one has to remember 
that pomeron does not change colour of colliding quarks) we have: 

4 ( A /e A ^) singlct 4 ( A fch A 4) singlet fsunfuwC = - — S fd 5 kg 5 uw , (18) 

what is a consequence of: 

2 

( A /e A ^) singlct = ^ AB5 fd 

WX 9 ) smglct = ls CA S k9 . (19) 

One can check that for the second diagram the colour factor is the same. 
The matrix element is given by (see Fig. 4): 
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Mj? = -y 2 5l f 5 uw j d*Q5 (qf) 5 (g 2 ) D (Q 2 ) D D (/ 2 ) 
u (h) 7^ /lil v u (p^ u (k 2 ) 7 A folvU (p 2 ) 
sVptf (h, -r u -q) —^~V P aX (q, ~r 2 , h) e\e% 

+V KS (h, ~r 2 , 0) ^Vpax (-9, ~ri, k) e?e|} , (20) 

where: 



^iiwiw ( fc i> fc 2, fc 3 ) = W2 (*i - fo)^ + g^ 3 (k 2 - h)^ + (k 3 - h)^ . (21) 

Here A = iD^G 4 g 2 /(2ir) 2 , G and g are the non-perturbative and perturbative quark- 
gluon couplings, D (p 2 ) = D exp (p 2 /fi 2 ) with \i ss 1 GeV and G 2 -D ~ 30 GeV _1 /i _1 
is the non-perturbative gluon propagator. The fact that we put inner quark lines on 
shell is expressed by 5 (q\) 5 (q 2 ). 

After some approximations, based on the fact that integrand in (20) is exponentially 
damped in v 2 , what allows to consider v as small, we obtain (for more details see [5]): 

Mf = -A^8%8 UW J d 2 vf (v) exp , 

where: 



/(#) = •{ p1V ks (h, -ri, -q) -^-V paX (q, -r 2 , h)p$e{e a 



+PiV^S (h, ~r 2 , q) (~q, ~n, Z 2 )z>2 e i£* 



Substituting (21) and (2) to the above formula, using p 2 = p\ = and S\T\ = e 2 r 2 
£iPi = ?2Pi = what is a consequence of (7), we obtain (6). 

Let us denote colour structure of the diagram represented by Afjf as follows: 



B i 



C i 




w 



Figure 8: The capital letters denote colour indicies of the non-perturbative gluons and 
small letters colour indicies of quarks and produced gluons. 

The matrix element Alif is given by (see Fig. 4): 
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< = ~4 I (^A^) singIet / d*Q6 (qf) 5 (q 2 2 ) 

D (Q 2 ) D D (l 2 2 ) u (ki) ^ kil v u (pi) u (k 2 ) 7 A fo>y v u (p 2 ) 

+ r Bu r Cw (9^a - gpagtx) } e\el. (22) 
Using (19) and performing the rest colour calculations we obtain: 

Mf = -^l f 6 uw j d*Q5 (qf) 5 (q 2 2 ) D (Q 2 ) D (Z?) D (l 2 2 ) 
u (h) I 13 h\l Vu (Pi) u (h) 7 A folvU (p 2 ) 

[29p\9z* - g^gxa - gp^x] 4 £ 2- ( 23 ) 

After some approximations, fully described in [5] we have: 



Mf = -A^5 i s { c S uw se 1 e 2 J d 2 vexp 
what is our result (4). 



-v 2 
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